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Abstract 

We present a measurement of the correlated hh production cross section. The data used in this 
analysis were taken with the upgraded CDF detector (CDF II) at the Fermilab Tevatron collider, 
and correspond to an integrated luminosity of 742 pb~^. We utilize muon pairs with invariant mass 
5 < m^j^^ < 80 GeV/c^ produced by hh double semileptonic decays. For muons with > 3 GeV/c 
and \r]\ < 0.7, that are produced by h and b quarks with pt > 2 GeV/c and \y\ < 1.3, we measure 
Cft^^ = 1549 lb 133 pb. We compare this result with theoretical predictions and previous 
measurements. We also report the measurement of ac^^^c^n, a by-product of the study of the 
background to bb production. 

PACS numbers: 14.65.Fy, 14.65.Dw, 12.38.Qk, 13.85.Qk, 13.20.He 
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I. INTRODUCTION 



Measurements of the cross section for producing, in hadronic collisions, both b and 
b quarks centrally and above a given transverse momentum threshold (typically pt > 
5—20 GeV/c) , referred to as (jj,^ or bb correlations, provide an important test of the predictive 
power of quantum chromo dynamics (QCD). Experimentally, bb correlations at the Tevatron 
are inferred from the production rate above a given px threshold of some of the decay prod- 
ucts (leptons or tracks consistent with a secondary displaced vertex) of both b and b hadrons. 
In QCD calculations, the long- and short-distance dynamics of the hadronic hard-scattering 
cross section are factorized into nonperturbative parton distribution functions (PDF) and 
fragmentation functions, and perturbatively calculable hard-scattering functions. At the 
perturbative level, the hard-scattering function can be evaluated at leading-order (LO) and 
next-to-leading order (NLO) with the mnr Monte Carlo program [1]. In contrast with the 
exact NLO prediction of the single b quark production cross section ^, the exact NLO cal- 
culation of ajji appears to be a robust perturbative QCD prediction. As noted in Ref. [5], 
the exact LO and NLO prediction of cr^^ are equal within a few percent, and the NLO re- 
sult does not change by more than 15% when varying the renormalization and factorization 
scales by a factor of two and the b quark pole mass (m^ = 4.75 GeV/c^) by 0.25 GeV/c^. 
The exact NLO prediction of cr^g is quite insensitive to the choice of PDF fits when they in- 
clude HERA data and yield a value of the QCD coupling strength consistent with LEP data 
{as{mz) — 0.118) [6-8]. However, when comparing to the data, the apparent robustness 
of the ajjjj calculation could be spoiled by the inclusion of nonperturbative fragmentation 
functions that connect 6-quark and 6-hadron distributions. Traditionally, data to theory 
comparisons use a fragmentation model based on the Peterson function [9] with the e pa- 
rameter set to 0.006 according to fits to e+e^ data [10]. However, as noted in Ref. [11], 
the Peterson fragmentation function has been tuned to the data in conjunction with LO 
parton-level cross sections evaluated with parton-shower event generators, and cannot be 
consistently convoluted with the exact NLO calculation. As an example, the FONLL cal- 

^ The single b quark cross section can be evaluated at exact NLO accuracy with the nde Monte Carlo gener- 
ator [2]. The calculation is affected by an uncertainty as large as 50% due to the choice of renormalization 
and factorization scales and by additional, but smaller, uncertainties due to choice of the PDF fits or the 
6-quark mass [3] . At perturbative level, the large scale dependence of the NLO calculation is interpreted 
as a symptom of large higher-order contributions [4] . 
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culation [12] implements the exact NLO prediction of the single 6-quark cross section with 
the resummation of {pr/iTib) logarithms with next-to-leading accuracy (NLL). A calculation 
with the same level of accuracy, available for the production of b quarks at e'^e~ colliders [13] , 
has been used to extract consistent nonperturbativc fragmentation functions from LEP and 
SLC data [14]. These fragmentation functions appear to be harder than the Peterson frag- 
mentation [11]. Unfortunately, they also cannot be consistently convoluted with the exact 
NLO calculation of cr^g for which NLL logarithmic corrections have yet to be evaluated. 

Alternatively, the production of pairs of b and b hadrons can be estimated with event 
generators that are based on the LO calculation combined with a leading-logarithmic (LL) 
treatment of higher orders via the parton shower approximation, such as the HERWIG [15] 
and PYTHIA [16] Monte Carlo programs. The MC@NLO event generator [17] merges the 
exact NLO matrix element with the LL shower evolution and hadronization performed by 
the HERWIG parton-shower Monte Carlo. In some cases, event generators that combine 
exact LO or NLO calculations with LL parton-shower simulations return parton-level cross 
sections that are quite different from the exact NLO calculation [1] . The mc@nlo method 
suffers the additional problem that the HERWIG model of the b quark hadronization has been 
tuned to e+e^ data using LO parton-level cross sections. The benefits and pitfalls of each 
theoretical approach are discussed in more detail in Refs. [17-19]. 

Precise measurements of the pair production of b and b hadrons at the Tevatron could 
contribute to improve the modeling of fragmentation functions consistent with the exact 
NLO calculation. Unfortunately, as noted in Ref . [6] , the status of the cr^g measurements at 
the Tevatron is quite disconcerting. Five measurements of af^i have been performed by the 
CDF and collaborations. Reference [6] compares the results of different experiments 
using i?2f>, the ratio of the measured aj,i to the exact NLO prediction (the 6-quark and b- 
hadron distributions are connected via the LL HERWIG fragmentation model or the Peterson 
fragmentation function) . 

The study in Ref. [5] (CDF) uses two central jets with Et > 15 GeV, each containing a 
secondary vertex due to b- or 6-quark decays. The measurement yields i?26 = 1-2 ± 0.3. 

The study in Ref. [20] (CDF) uses events containing two central jets with Et > 30 and 
20 GeV, respectively; pairs of b jets are also identified by requiring the presence of displaced 
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secondary vertices. This study yields ^ i?26 = 1.1 ± 0.3. 

The study in Ref. [21] (CDF) uses events containing muons from 6-quark semileptonic 
decays that recoil against a jet that contains tracks with large impact parameter {b jet). 
This study yields i?26 = 1.5 ±0.2 for b and b quarks produced centrally with > 12 GeV/c. 

References [22] (CDF) and [24] (DjZ)) report measurements that use two central muons 
arising from 6-quark semileptonic decays. The measurements yield i?26 = 3.0 ± 0.6 and 
i?26 = 2.3 ± 0.7 for central b and b quarks with pt > 6 and 7 GeV/c, respectively. 

The five measurements yield < i?26 >= 1-8 with a 0.8 RMS deviation [6]. Such a large 
RMS deviation is a likely indication of experimental difficulties ^. This type of discrepancy 
could result from an underestimate of the kinematic and detector acceptance for semileptonic 
b decays or of the underlying background. However, measurements of the single 6-quark pro- 
duction cross section based upon detection of semileptonic 6-quark decays suggest otherwise 
because they are approximately 35% smaller than those based on detection of J/ip mesons 
from 6-quark decays [6, 23] . The present discrepancy could also be explained by postulating 
the production of additional objects with a 100% semileptonic branching ratio and a cross 
section of the order of 1/10 of the b cross section as investigated in Ref. [5]. Therefore, it 
is of interest to clarify the experimental situation. This paper reports a new measurement 
of (7^5 that uses dimuons arising from bb production. At the Tevatron, dimuon events result 
from decays of heavy quark pairs {bb and cc), the Drell-Yan process, charmonium and bot- 
tomonium decays, and decays of tt and K mesons. Background to dimuon events also comes 
from the misidentification of tt or X mesons. As in previous studies [22, 25], we make use 
of the precision tracking provided by the CDF silicon microvertex detector to evaluate the 
fractions of muons due to long-lived b- and c-hadron decays, and to the other background 
contributions. 

Sections II and III describe the detector systems relevant to this analysis and the data 
selection, respectively. The analysis method is discussed in Sec. IV, while the heavy flavor 
composition of the dimuon sample is determined in Sec. V. The kinematic and detector 
acceptance is evaluated in Sec. VI. The dimuon cross section is derived and compared 
to theoretical expectation and previous measurements in Sec. Vll. Our conclusions are 

^ Ref. [20] compares the data to the mc@nlo prediction, which is 12% smaller than the exact NLO predic- 
tion. This difference was not appreciated in Ref. [6]. 
^ This includes the possibility that in some cases the NLO prediction has been evaluated incorrectly. 
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summarized in Sec. VIII. 



II. CDF II DETECTOR AND TRIGGER 

CDF II is a multipurpose detector, equipped with a charged particle spectrometer and a 
finely segmented calorimeter. In this section, we describe the detector components that are 
relevant to this analysis. The description of these subsystems can be found in Refs. [26-35] . 
Two devices inside the 1.4 T solenoid are used for measuring the momentum of charged 
particles: the sihcon vertex detector (SVXII and ISL) and the central tracking chamber 
(COT). The SVXII detector consists of microstrip sensors arranged in six cyhndrical shells 
with radii between 1.5 and 10.6 cm, and with a total z coverage ^ of 90 cm. The first 
SVXII layer, also referred to as LOO detector, is made of single-sided sensors mounted on the 
beryllium beam pipe. The remaining five SVXII layers are made of double-sided sensors and 
are divided into three contiguous five- layer sections along the beam direction z. The vertex 
2;-distribution for pp collisions is approximately described by a Gaussian function with a 
sigma of 28 cm. The transverse profile of the Tevatron beam is circular and has an RMS 
spread of ~ 25 \xm. in the horizontal and vertical directions. The SVXII single-hit resolution 
is approximately 11 /xm and allows a track impact parameter resolution of approximately 
35 /xm, when also including the effect of the beam transverse size. The two additional 
sihcon layers of the ISL help to link tracks in the COT to hits in the SVXII. The COT is 
a cylindrical drift chamber containing 96 sense wire layers grouped into eight alternating 
superlayers of axial and stereo wires. Its active volume covers \z\ < 155 cm and 40 to 140 
cm in radius. The transverse momentum resolution of tracks reconstructed using COT hits 
is ct{pt)/Pt — 0.0017 [GcV/c]^^. COT tracks are extrapolated into the SVXII detector and 
refitted adding hits consistent with the track extrapolation. 

The central muon detector (CMU) is located around the central electromagnetic and 
hadronic calorimeters, which have a thickness of 5.5 interaction lengths at normal incidence. 

In the CDF coordinate system, 9 and are the polar and azimuthal angles of a track, respectively, defined 
with respect to the proton beam direction, z. The pseudorapidity rj is defined as — log ta.n{9/2). The 
transverse momentum of a particle is pt = p sin(^). The rapidity is defined as y = l/2-log{{E + pz)/{E — 
Pz)), where E and Pz are the energy and longitudinal momentum of the particle associated with the track. 
^ The impact parameter d is the distance of closest approach of a track to the primary event vertex in the 
transverse plane. 
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The CMU detector covers a nominal pseudorapidity range I77I < 0.63 relative to the center of 
the detector, and is segmented into two barrels of 24 modules, each covering 15° in 0. Every 
module is further segmented into three submodules, each covering 4.2° in cf) and consisting 
of four layers of drift chambers. The smallest drift unit, called a stack, covers a 1.2° angle 
in (p. Adjacent pairs of stacks are combined together into a tower. A track segment (hits 
in two out of four layers of a stack) detected in a tower is referred to as a CMU stub. A 
second set of muon drift chambers (CMP) is located behind an additional steel absorber of 
3.3 interaction lengths. The chambers are 640 cm long and are arranged axially to form a 
box around the central detector. The CMP detector covers a nominal pseudorapidity range 
\ri\ < 0.54 relative to the center of the detector. Muons which produce a stub in both CMU 
and CMP systems are called CMUP muons. 

The luminosity is measured using gaseous Cherenkov counters (CLC) that monitor the 
rate of inelastic pp coUisions. The inelastic pp cross section at ^Js — 1960 GeV is scaled 
from measurements at ^/s — 1800 GeV using the calculations in Ref. [36]. The integrated 
luminosity is determined with a 6% systematic uncertainty [37]. 

CDF uses a three-level trigger system. At Level 1 (LI), data from every beam crossing 
are stored in a pipeline capable of buffering data from 42 beam crossings. The LI trigger 
either rejects events or copies them into one of the four Level 2 (L2) buffers. Events that pass 
the LI and L2 selection criteria are sent to the Level 3 (L3) trigger, a cluster of computers 
running speed-optimized reconstruction code. 

For this study, we select events with two muon candidates identified by the LI and L2 
triggers. The LI trigger uses tracks with pr > 1.5 GeV/c found by a fast track processor 
(XFT). The XFT examines COT hits from the four axial superlayers and provides r — 
information. The XFT finds tracks with px > 1.5 GeV/c in azimuthal sections of 1.25°. 
The XFT passes the tracks to a set of extrapolation units that determine the CMU towers 
in which a CMU stub should be found if the track is a muon. If a stub is found, a LI 
CMU primitive is generated. The LI dimuon trigger requires at least two CMU primitives, 
separated by at least two CMU towers. The L2 trigger additionally requires that at least 
one of the muons has a CMUP stub matched to an XFT track with > 3 GeV/c. All these 
trigger requirements are emulated by the detector simulation on a run-by-run basis. The 
L3 trigger requires a pair of CMUP muons with invariant mass larger than 5 GeV/c^, and 
\5zo\ < 5 cm, where zq is the z coordinate of the muon track at its point of closest approach 
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to the beam line in the r — plane. These requirements define the dimuon trigger used in 
this analysis. 

We use additional triggers in order to measure detection efficiencies and verify the detector 
simulation. The first trigger (CMUPpr4) selects events with at least one CMUP primitive 
with Pt > 4: GeV/c identified by both the LI and L2 triggers, and an additional muon 
found by the L3 algorithms. Events collected with this trigger are used to measure the muon 
trigger efficiency. The second trigger requires a LI CMUP primitive with pt > 4 GeV/c 
accompanied by a L2 requirement of an additional track with Pt > 2 GeV/c and impact 
parameter 0.12 < c? < 1 mm as measured by the Sihcon Vertex Trigger (SVT) [38]. The 
SVT calculates the impact parameter of each XFT track, with respect to the beam line, 
with a 50 fira resolution that includes the 25 /xm contribution of the beam transverse width. 
Events selected with this trigger (/x— SVT) are used to verify the muon detector acceptance 
and the muon reconstruction efficiency The last trigger (charm) acquires events with two 
SVT tracks with pt > 2 GeV/c and with impact parameter 0.12 < d < 1 mm. In this data 
sample, we reconstruct D'^ — > Ktt decays to measure the probability that a charged hadron 
mimics the signal of a CMUP muon. We also use J/ip —>■ fi'^fi~ events acquired with the 
J/iJj trigger. At LI and L2, this trigger requires two CMU primitives corresponding to tracks 
with pt > 1.5 GeV/c. At L3, muons are required to have opposite charges and an invariant 
mass in the window 2.7 — 4.0 GeV/c^. These events are used to calibrate the efficiency of 
the SVXII detector and of stricter requirements used for selecting CMUP muons. 

III. DATA SELECTION 

In this analysis, we select events acquired with the dimuon trigger and which contain two 
and only two CMUP muons with same or opposite charge. Events are reconstructed offline 
taking advantage of more refined calibration constants and reconstruction algorithms. COT 
tracks are extrapolated into the SVXII detector, and refitted adding hits consistent with 
the track extrapolation. Stubs reconstructed in the CMU and CMP detectors are matched 
to tracks with Pt > 3 GeV/c. A track is identified as a CMUP muon if Ar0, the distance 
in the r — plane between the track projected to the CMU (CMP) chambers and a CMU 
(CMP) stub, is less than 20 (40) cm. We require that muon- candidate stubs correspond to 
a LI CMU primitive, and correct the muon momentum for energy losses in the detector. 
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To ensure an accurate impact parameter measurement, each muon track is required to 
be reconstructed in the SVXII detector with hits in the two inner layers and in at least 
two of the remaining four external layers. We evaluate the impact parameter of each muon 
track with respect to the primary vertex. We reconstruct primary vertices using all tracks 
with SVXII hits that are consistent with originating from a common vertex. In events in 
which more than one interaction vertex has been reconstructed we use the one closest in z to 
the average of the muon track 2;o-positions and within a 6 cm distance. The primary vertex 
coordinates transverse to the beam direction have RMS uncertainties of approximately 3 //m, 
depending on the number of SVXII tracks associated with the primary vertex and the event 
topology. 

Muon pairs arising from cascade decays of a single h quark are removed by selecting 
dimuon candidates with invariant mass greater than 5 GeV/c^. We also reject muon pairs 
with invariant mass larger than 80 GeV/c^ that are mostly contributed by decays. 

IV. METHOD OF ANALYSIS 

For muons originating from the decay of long lived particles, the impact parameter is 
d = \f3^ct sm{6)\, where t is the proper decay time of the parent particle from which the 
muon track originates, 6 is the decay angle of the muon track with respect to the direction of 
the parent particle, and is the Lorentz boost factor. The impact parameter distribution 
of muon tracks is proportional to the lifetime of the parent particle. The markedly different 
distributions for muons from b decays, c decays, and other sources allow the determination 
of the parent fractions. 

We determine the bb and cc content of the data following the method already used in 
Refs. [22, 25]. The procedure is to fit the observed impact parameter distribution of the muon 
pairs with the expected impact parameter distributions of leptons from various sources. After 
data selection, the main sources of reconstructed muons are semileptonic decays of bottom 
and charmed hadrons, prompt decays of quarkonia, and Drell-Yan production. 

Monte Carlo simulations are used to model the impact parameter distributions of muons 
from b- and c-hadron decays. We use the HERWIG Monte Carlo program [15], the settings 
of which arc described in Appendix A, to generate hadrons with heavy flavors that are 
subsequently decayed using the evtgen Monte Carlo program [39]. The detector response 
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to particles produced by the above generators is modeled with the CDF II detector simu- 
lation that in turn is based on the GEANT Monte Carlo program [40]. Impact parameter 
distributions of muon tracks in simulated b- and c-hadron decays are shown in Fig. 1. Since 
lifetimes of bottom and charmed hadrons (cr^ ~ 476 and ctq — 213 /im) are much 
larger than the average SVXII impact parameter resolution (~ 28 /im), the dominant factor 
determining the impact parameter distribution is the kinematics of the semileptonic de- 
cays which is well modeled by the evtgen program. The impact parameter distribution 
of muons from prompt sources, such as quarkonia decays and Drell-Yan production, is con- 
structed using muons from T(IS') decays (see Fig. 2). Muons from tt and K in-flight decays 
are also regarded as prompt tracks since the track reconstruction algorithm rejects those 
with appreciable kinks. Tracks associated with tt and K mesons which mimic the lepton 
signal (fake muons) are mostly prompt. The small contribution to fake muons of pion and 
kaon tracks arising from the decay of hadrons with heavy flavor is evaluated separately in 
Sec. V B. Since there are two muons in an event, the flt is performed in the two-dimensional 
space of impact parameters. Each axis represents the impact parameter of one of the two 
muons. In filling the histograms, the muon assignment is randomized. The two-dimensional 
impact parameter technique exploits the fact that muon impact parameters are independent 
uncorrelated variables ^. The two-dimensional template distributions for each type of event 
are made by combining the relevant one-dimensional distributions in Fig. 1. 

We use a binned maximum log likelihood method [41] to fit the dimuon impact parameter 
distribution. The hkelihood function L is defined as 

L-iimf'^^-'vri{z,m (1) 

i j 

where n{i,j) is the number of events in the (i, j)-th bin. The function l^j is defined as 

kj = BB ■ S,{i) ■ S,{j) + CC ■ S,{i) ■ S,{j) + PP ■ S,{i) ■ S,{j) + (2) 
0.5 • [BP ■ {S,{{} ■ Sp{j) + Spi{} ■ S,{j)) + CP ■ {S,{{) ■ Sp{j) + Sp{{) ■ S,{j)) + 
BC ■ {S,{t) ■ S,{j) + S,{t) ■ S,{j))] 

where Sb, Sc, and Sp are the impact parameter templates shown in Fig. 1. The fit parameters 
BB, CC, and PP represent the bb, cc, and prompt dimuon contributions, respectively. The 



^ The correlation between the two impact parameters, p = ^ — ^ ^^J^ 1 ^ — ^- — -) is approximately 
0.01 in the data and the heavy flavor simulation. 
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FIG. 1: Impact parameter distributions of muons coming from b- and c-hadron decays (simulation) 
and of prompt muons (data) . Distributions are normalized to unit area. 



fit parameter BP (CP) estimates the number of events in which there is only one b (c) 
quark in the detector acceptance and the second lepton is produced by the decay or the 
misidentification of vr or mesons The fit parameter BC estimates the number of events 
in which both bottom and charmed quarks are final state partons of the hard scattering. 

^ According to the simulation, approximately 86% of the bb and cc events with an identified muon from 
heavy flavor decay do not contain a second hadron with heavy flavor in the detector acceptance. Therefore, 
following the procedure of Ref. [25], we start by ignoring the small fake muon contribution due to tt and 
K mesons from heavy flavor decays, that is estimated in Sec. VB. 
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FIG. 2: Distribution of the invariant mass of muon pairs in the T region. The prompt template in 
Fig. 1 is derived using muons with invariant mass between 9.28 and 9.6 GeV/c^. The background is 
sideband subtracted using dimuons with invariant mass between 9.04 and 9.2 GeV/c^ and between 
9.64 and 9.8 GeV/c^. 

According to the simulation, the BC component is ~ 4.6% of the BB component and the 
CP component is ~ 83% of the BP component Figure 3 shows projections of the two- 

® In the simulation, events containing a muon from heavy flavor decay and a prompt track are mostly 
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dimensional distributions for each type of mixed contribution. By comparing with Fig. 1, 
one notes that the BB and PP components have impact parameter distributions markedly 
different from any other contribution, whereas the CC, CP, BP, and BC components 
have quite similar shapes. Using Monte Carlo pseudoexperiments we have verified that, as 
observed in previous studies [22, 25], the likelihood function is not capable of disentangling 
these four components. Therefore, Eq. (1) is supplemented with the term 

(CP - 0.83 -^P)^ {BC-0M6-BBy 
■ ' ^ {CP + 0.832 . BP + (0.14 • BPy) {BC + 0.046^ • BB + (0.013 • BB)^) ' ^ ' 

that constrains the ratios CP/ BP and BC/BB to the values predicted by the simulation 
within their theoretical uncertainties approximated with Gaussian functions ^. 
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FIG. 3: Projections of the two-dimensional impact parameter distributions of some components 
used to fit the dimuon data (see text). 



contributed by NLO diagrams, such as flavor excitation and gluon splitting, in which a heavy flavor quark 

recoils against a gluon or a light quark. The cross section for producing at least one c quark in the 
kinematic acceptance of this study is 2.6 times larger than that of a & quark, but the contribution of NLO 
terms in cc production is approximately 3.6 times larger than for bb production. However, the kinematic 
acceptance for muons from c decays is ~ 23% of that for b decays. 

Using other PDF fits available in the PDF library [42], the c-to-b ratio of the flavor excitation cross section 

in the simulation varies up to ±30% [43]. The ratio of the c-to-b gluon splitting cross section also changes 
by ±30% when varying the c- and 6-quark pole mass by 0.5 GcV/c^ [43]. We use as Gaussian uncertainty 
the 60% variation divided by The ratio of be to bb production depends on the c-quark structure 

function, and varies up to ±50% when using other PDF fits [42]. 
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V. HEAVY FLAVOR COMPOSITION OF THE DIMUON SAMPLE 



In this section, we first determine the dimuon sample composition by fitting the impact 
parameter distribution with the templates described in the previous section. We then eval- 
uate and remove the contribution of muons faked by pion or kaon tracks from heavy flavor 
decays. Lastly, we estimate the systematic uncertainty of the result due to the fit likelihood 
function and simulated templates. 



A. Result of the fit to the impact parameter distribution 

The two-dimensional impact parameter distribution of the 161948 muon pairs selected 
in this analysis is plotted in Fig. 4(a). An appreciable fraction of events cluster along the 




FIG. 4: Two-dimensional impact parameter distributions of muon pairs (a) before and (b) after 
cosmic removal. 



diagonal line di = d2- These events are due to cosmic rays and wc remove them by requiring 
the azimuthal angle between muons with opposite charge to be smaller than 3.135 radians 
(see Fig. 4(b)). In the simulation, the Scf) < 3.135 requirement has a 99.3% efficiency when 
applied to muon pairs arising from bb and cc production. The invariant mass spectrum of 
the remaining 143677 events is shown in Fig. 5. 
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FIG. 5: Invariant mass spectrum of the muon pairs used in this study. 



The result of the fit to the two-dimensional impact parameter distribution of the data 
using the likehhood function in Eqs. (1-3) is shown in Table I. The parameter correlation 
matrix is listed in Table II. The projection of the two-dimensional impact parameter 
distribution is compared to the fit result in Fig. 6 and the distribution of the fit residuals 

is plotted in Fig. 7. The best fit returns — InL = 1078. The probability of the fit to the 
data is determined with Monte Carlo pseudoexperiments. In each experiment, we randomly 
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TABLE I: Number of events attributed to the different dimuon sources by the fit to the impact 
parameter distribution. The errors correspond to a 0.5 change of — InL. 



Component No. of Events 

BB 54583 ± 678 

CC 24458 ± 1565 

PP 41556 ± 651 

BP 10598 ± 744 

CP 10024 ± 1308 

BC 2165 ± 693 



TABLE II: Parameter correlation coefficients returned by the fit hsted in Table I. 

Component BB CC PP BP CP 

CC -0.46 

PP 0.09 0.18 

BP 0.01 -0.43 -0.14 

CP 0.27 -0.69 -0.71 0.13 

BC -0.42 -0.19 0.15 -0.18 -0.06 



generate different components with average size as determined by the fit to the data 
and allow for Poisson fluctuations; the impact parameter distribution for each component 
is randomly generated from the corresponding templates used to fit the data. We find that 
16.5% of the fits to the pseudoexperiments return a — InL value equal or larger than 1078. 
The values of the different components returned by the fits to each pseudoexperiment have 
Gaussian distributions with sigmas equal to the corresponding errors listed in Table I. The 
fit result is quite insensitive to the constraint CP/BB — 0.83 ± 0.14 in Eq. (3). If the 
uncertainty is increased from 0.14 to 0.28, the size of the BB and CC components returned 

We have performed 1000 pseudoexperiments starting from the following component sizes: BBq = 54600, 
CCo = 24500, PPo = 41500, BPq = 10200, CPq = 10000, and BCo = 2200. 
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by the fit changes by less than half of the corresponding errors listed in Table I. However, 
without this constraint, the fit returns a CC component 30% smaller than the standard fit 
together with a ratio CP/ BP ~ 3.5 that would be difficult to account for. In comparison 




FIG. 6: The projection of the two-dimensional impact parameter distribution of muon pairs onto 
one of the two axes is compared to the fit result (histogram). 



with a previous CDF measurement [25] that uses data collected in the 1992 — 1995 collider 
run (Run I), the ratio PP/BB returned by the fit has increased from 34% to 76%. In the 
present study, we do not remove T candidates (13800 ±290) which represent 32% of the PP 
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FIG. 7: Distribution of the pulls - (data-fit)/-\/fit - of the fit listed in Table I for the impact 
parameter bins with at least 10 entries. The solid line represents a Gaussian distribution with unit 
RMS. 

contribution. After removing the T contribution, the PP contribution in the present data 

sample is 50% larger than in the Run I data. This is explained by a substantial increase in 
the rate of muons faked by prompt tracks. Therefore, in the next subsection, we evaluate 
the fraction of the BB yield due to muons faked by hadronic tracks from heavy flavor decays 
(this contribution was found negligible in the Run I data [22, 25]). 
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B. Fake muon contribution 



We use two methods to estimate the contribution of tracks arising from heavy flavor 
decays that mimic a CMUP signal. The first method is based on a combination of data and 
simulation. We use the simulation to estimate the relative yields, Rk and Rt^, ol ji — K 
and n — 7T combinations with respect to that of real muon pairs in the decay of hadrons 
with heavy flavor (we select muons and tracks with pr > 3 GeV/c and \ri\ < 0.7, and we 
require the invariant mass of each pair to be larger than 5 GeV/c^). These yields are listed 
in Table III. The ratio of dimuons contributed by //—track combinations to real dimuons 
from semileptonic decays is F = {Rk ■ Pf + -Rtt • Pj)/e^, where Pf^ (PJ) is the probability 
that a pion (kaon) track mimics a muon signal. These probabilities are determined using 
a sample of nK decays in App. B. The corresponding efficiency for detecting a real 

muon is = 0.5057 (see Table VIII). The errors in Tabic III arc the sum in quadrature of 
statistical errors and the 10% systematic uncertainty of the kaon and pion rates predicted 
by the simulation. This systematic uncertainty is derived from a comparison of kaon and 
pion production rates measured at the T(45') to the prediction of the evtgen Monte Carlo 
program [44]. 

TABLE III: Ratio of the numbers of /x — K{Tr) combinations to that of /x — // pairs, Rk{tt)-, in the 
simulation of different heavy flavor productions. The ratio F of the number of fake-real muon pairs 
to that of real dimuons is estimated using the fake muon probabihties derived in App. B and the 
measured detector efficiency = 0.5057 for a real muon. 



Production 

bb 

cc 

be 



Rk 
3.70 ± 0.43 
21.73 ±2.51 
16.78 ± 2.71 



Rw 
7.58 ± 0.82 
23.47 ± 2.68 
10.83 ± 2.03 



F (%) 
7.2 ±0.6 

32.4 ±2.7 

21.5 ±2.8 



We evaluate the purity (1/(1 + F)) with a second method that is almost independent of 
the simulation prediction. We make use of stricter muon selection criteria by supplementing 
the Ar(f) cut between the muon track projection and the CMU and CMP stubs with the 
requirement, referred to as cut, that the extrapolated COT track and the CMU muon 
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FIG. 8: Invariant mass distribution of CMUP muon pairs with a randomly chosen muon that (a) 
satisfies or (b) fails the < 9 requirement. 



stub match within 3 a in the r — plane, where a is a standard deviation that includes the 
effect of multiple scattering and energy loss. The efficiency of the cut for real muons is 
measured using a sample of muons acquired with the J/ ip trigger. We compare the invariant 
mass distributions of CMUP pairs when a randomly chosen muon passes or fails the ^ 9 
cut. We fit the data with two Gaussian functions to model the J/ip signal and a straight 
line to model the background. The rj and pt distribution of CMUP muons from J/ ip decays 
are weighted to model that of muons from 6-hadron decays. As shown in Fig. 8, the cut 
reduces the efficiency for detecting a muon pair by einefr = 2.20 ± 0.04%. 

The corresponding fake muon probabilities are measured using Ktc decays and 

are listed in Table XI of App. B. We select a sample of dimuons enriched in fake muons by 
requiring > 9 for one muon and determine its heavy flavor composition by fitting the 
impact parameter distribution. The fit result is shown in Table IV. For each heavy fiavor 
component, we derive the fake muon contribution by solving the system of equations 

T = HF + PfFK (4) 
T{x' > 9) = eineff ■ HF + Pf{x' > 9) ■ FK, 
where T and T(x^ > 9) are the size of the component determined by the fits in Table I 
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TABLE IV: Number of events attributed to the different dimuon sources by the fit to the impact 
parameter distribution. We use events in which at least one muon fails the > 9 requirement 
(see text). 



Component No. of Events 

BB 1103 ± 102 

CC 1189 ± 272 

PP 4249 ± 131 

BP 1508 ± 136 

CP 1218 ± 194 

BC 51 ± 15 



and IV, respectively, HF is the number of real muon pairs, and FK is the number of 
dimuons one of which is faked by a track from heavy flavor decays. The fraction of real 
muon pairs reads 

^ ' T • {Pfix^ > 9) - 0.022 • Pf) 



± r . (PKx- > 9) - 0.022 . p,) > °) + ^'-^^ > ■ 

This second method provides a determination of the fraction of real dimuons almost inde- 
pendent of the pion and kaon rate predicted by the simulation. The fraction of real muon 
pairs determined with the two methods is shown in Table V. We use the average and take 
the maximum and minimum RMS deviation as systematic uncertainty (0.96 ± 0.04 for bb 
and 0.81 ± 0.09 for cc production). The contribution of pairs of muons that are both faked 
by tracks from heavy flavor decays has been estimated to be less than 0.4% in the worst 
case, and it is ignored. 

C. Results after fake removal 

Table VI lists the various heavy flavor contributions to the dimuon sample after removing 
the contribution of tracks from heavy flavor decays that mimic a muon signal. As shown in 
Table V, the contribution of muons faked by tracks from c-quark decays is not negligible. 
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TABLE V: Fractions of real dimuons due to heavy flavor, 1/(1+F), determined with the simulation 
or by using the results returned by fits to the impact parameter distributions of all muon pairs and 
of those pairs in which at least one muon fails the < 9 cut. 



Production 
bb 



Simulation 
0.93 ± 0.01 
0.76 ± 0.02 



>9 
1.01 ± 0.01 
0.86 ± 0.06 



Therefore, we search the simulation for combinations of muons from b semileptonic decays 
and pion or kaon tracks from b- or c-quark decays (both with Pt > 3 GeV/c and \r]\ < 0.7). 
Figure 9 compares impact parameter distributions of these tracks to the standard muon 
templates. Distributions for muons and hadrons are quite similar. We have fitted the data 
with templates that include the expected contribution of muons faked by tracks from heavy 
flavor decays as listed in Table V. The result of this fit differs by less than 0.1% from that 
of the standard fit in Table I. 

TABLE VI: Number of real muon pairs from heavy flavor sources after removing the fake muon 
contributions. Errors include the uncertainty of the fake removal. 



Component 
BB 

CC 



No. of Events 
52400 ± 2278 
19811 ± 2540 



For completeness, we use the simulation to verify the ratio of the BP to BB components 
returned by the fit performed in Sec. V. The fit yields a ratio BP/BB = 0.194 ± 0.013. 
We search the simulation for combinations of muons from b semileptonic decays and prompt 
pion or kaon tracks (both with pt > 3 GeV/c and |?7| < 0.7). The ratio of their number 
to that of dimuons from b semileptonic decays is 32.8 ± 0.6 (stat.). Since the efficiency for 
detecting a muon is = 0.5 and the probability that prompt tracks fake a muon signal is 
0.0032, the simulation predicts BP/BB — 0.21 ± 0.01, in fair agreement with the fit result 
even without considering the uncertainty of the rate of prompt pions and kaons predicted 
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d (cm) d (cm) 

FIG. 9: Simulated impact parameter distributions of pion and kaon tracks from (a) b- and (b) c- 
quark decays are compared to those of muons from semileptonic decays with the same kinematical 
requirements. 

by the simulation. In cc data, approximately 20% of the muons are faked by hadronic tracks 
from c-quark decays. This is in agreement with the result that CP/ BP is 0.95 ± 0.14 in the 
data and 0.83 in the heavy flavor simulation. 

D. Dependence of the result on the muon px distribution 

The impact parameter of a track arising from heavy flavor decays depends on the proper 
decay time of the parent hadron and the decay angle between the daughter track and the 
parent hadron in its rest frame. Because we select muons above a given pt threshold, the 
range of accepted decay angles shrinks as the pr difference between the daughter track and 
parent hadron decreases. Therefore, impact parameter templates have a small dependence on 
the transverse momentum distribution of muons (or parent heavy flavor) in the simulation. 

In the data, we derive transverse momentum distributions for muons from h- and c-hadron 
decays by using the sP/ot statistical method [45] . We call /„ one of the six components used 
in the likelihood function L in Eq. (2), Nn the number of events attributed to this component 
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FIG. 10: Transverse momentum distributions in the data (•) and simulation (o) for muon pairs 
arising from (left) bb and (right) cc production. Distributions are normalized to unit area. 

by the fit in Table I, and the total number of events Given an event e in which muons 
have impact parameters in the (ze,Je)-th bin, the probability that the event belongs to the 
n-th component is 



(6) 



where 



Tr — 1 X]e=l /n(^e; Je) ' fli^ei je) 



(7) 



The transverse momentum distribution of muons from bb and cc production is obtained by 
weighting the muon transverse momenta in the event e by the corresponding probabilities 
Pi{ie,je) and P2{ie,je), respectively. The corresponding errors have been evaluated with 
Monte Carlo pseudoexperiments. These distributions are compared to those in the simula- 
tion in Fig. 10. The bb data are quite well modeled by the HERWIG generator. We estimate 
the dependence of the fit result on the muon px spectrum in the heavy flavor simulation by 
fitting the function A ■ to the ratio of these distributions in the simulation and in the 
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For example, fi = Sb- Sb and A^i = BB. 
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FIG. 11: Transverse momentum distribution of (left) prompt muons in the data and (right) muons 
from b- or c-hadron decays in the simulation. Simulated distributions are normalized to unit area. 



data. The fit returns a = — 0.029±0.015. We have constructed templates by reweighting the 
simulated muon transverse momentum distribution with the function p^^ *^^^. When these 
templates are used, the BB and CC yields returned by the fit change by ±1.5%, and =f4%, 
respectively. We do not correct our result for this effect, but we add this variation to other 
systematic effects evaluated in Sec. V E. 

For completeness, we use the probability defined in Eq. (6) to show: a comparison of 
transverse momentum distributions of prompt muons and muons from simulated heavy flavor 
decays in Fig. 11; distributions of 6(j), the azimuthal opening angle between two muons, in 
the data and the simulation in Fig. 12; the invariant mass spectrum of the PP component 
in the T mass region in Fig. 13; and a data to simulation comparison of the invariant mass 
spectrum of dimuon pairs from heavy flavor decays in Fig. 14. 



E. Dependence of the result on the b- and c-quark lifetime 

The lifetime of the b hadron mixture with semileptonic decays produced at the Tevatron 
has a 0.6% uncertainty [47]. Impact parameter templates, constructed by varying the lifetime 
by this uncertainty, change the BB size returned by the fit by ±0.4% and the CC size by 
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FIG. 12: Distributions of the opening azimuthal angle between two muons from (left) bb and (right) 
cc production in the data (•) and simulation (o). Distributions are normalized to unit area. 






9 9.5 10 10.5 

m^^ ( GeV/c' ) 

FIG. 13: Distribution of the invariant mass of the prompt dimuon component as determined by 
the fit to the muon impact parameter. The solid line represents a fit to the distribution that 
returns 9899 it 142 T{1S) candidates, whereas the same fit to the data in Fig. 2 yields 9952 it 122 
candidates. 
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FIG. 14: Distributions of the invariant mass of muon pairs from (left) bb and (right) cc production 
in the data (•) and simulation (o). Distributions are normalized to unit area. 



±1%. The Ufetime of the c-hadron mixture has a ±3.2% uncertainty, mostly due to the 
uncertainty of the relative fractions of produced hadrons, listed in Table X, and to the 
uncertainty of the semileptonic branching fractions of different c hadrons [47]. When using 
simulated templates constructed by changing the average lifetime by ±3.2%, the BB size 
returned by the fit changes by ±1% and the CC size varies by ±3%. By adding linearly 
these systematic uncertainties to that due to the muon px spectrum in the simulation, we 
derive a ±2.9% systematic error for the BB component and ±8% systematic error for the 
CC component. 



VI. ACCEPTANCE AND EFFICIENCIES 



The kinematic and detector acceptance is calculated with the Monte Carlo simulation 
described at the beginning of Sec. IV and in App. A. The detector response to muons 
produced by b- and c-hadron decays is modeled with the CDF II detector simulation that 
also models the LI and L2 trigger responses. Simulated events are processed and selected 
with the same analysis code used for the data. The acceptance (A) is the fraction of 
generated muon pairs that are identified in the detector and pass all selection requirements. 
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TABLE VII: Detector and kinematic acceptances, A, for dimuon pairs arising from bb and cc 
production. The acceptance .4corr includes corrections evaluated using the data. 



Production 
bb 



A (%) 
4.21 ± 0.04 
3.95 ± 0.10 



"^corr (%) 

4.56 ±0.15 
4.28 ±0.17 



At generator level, we select pairs of muons with invariant mass 5 < runn < 80 GeV/c^, each 
having Pt > 3 GeV/c and I77I < 0.7. Acceptances derived from the simulation are listed 
in Table VII. We use the data to verify the detector acceptance and efficiencies evaluated 
using the CDF II detector simulation. We adjust the simulation to match measurements 
in the data of: (1) the offline COT track reconstruction efficiency; (2) the CMUP detector 
acceptance and efficiency; (3) the efficiency for finding LI CMU primitives; (4) the SVXII 
acceptance and efficiency; and (5) the efficiency of the LI, L2, and L3 triggers. 

In the simulation, the offiine COT track reconstruction efficiency (0.998 ± 0.002) is the 
fraction of tracks, which at generator level satisfy the pr and rj selection cuts, that survives 
after selecting fully simulated events as the data. In the data, this efficiency has been 
measured to be 0.996 with a ~ 0.006 systematic accuracy by embedding COT hits generated 
from simulated tracks into J/i/j data [58] As in a previous study [23], we conclude that 
the efficiencies for reconstructing muon pairs in the data and the simulation are equal within 
a 1.3% systematic uncertainty. 

In the simulation, the fraction of CMUP stubs generated by muon tracks with pt > 
3 GeV/c and \ri\ < 0.7 is 0.5235 ± 0.0022. In the data, this efficiency is measured by using 
J/ip fi~^ fi~ decays acquired with the ^-SVT trigger. We evaluate the invariant mass of 
all pairs of a CMUP track and a track with displaced impact parameter, pr > 3 GeV/c, 
and I77I < 0.7. We fit the invariant mass distribution with a first order polynomial plus two 
Gaussian functions to extract the J/ip signal. Prom the number of J/ip mesons reconstructed 
using displaced tracks with or without a CMUP stub (Fig. 15(a) and (b), respectively), we 

The efficiency measurement was performed in a subset of the data used for this analysis. Studies of 
independent data samples collected in the data taking period used for this analysis show that changes of 
the track reconstruction efficiency are appreciably smaller than the quoted systematic uncertainty [59]. 
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FIG. 15: Invariant mass distribution of CMUP muons paired with a displaced track with (a) or 
without (b) a CMUP stub. Lines represent the fits described in the text. 



derive an efficiency of 0.5057 ± 0.0032. The integrated efficiency is evaluated after having 
weighted the px and r] distributions of displaced tracks in the data to be equal to those of 
muons from heavy flavor decays in the simulation. 

In the simulation, the efficiency for finding a LI CMU primitive (CMU stub matched by 
a XFT track) is 0.8489 ± 0.0026. This efficiency is measured in the data by using events 
acquired with the CMUPpt4 trigger. We combine the CMUP muon with all other CMUP 
muons found in the event with and without a LI CMU primitive. We extract the number 
oi J/ip —* mesons by fitting the invariant mass distributions of all candidates with 

a first order polynomial plus two Gaussian functions. By comparing the fitted numbers of 
J/ip candidates with and without a LI CMU primitive (Fig. 16(a) and (b), respectively) we 
derive an efficiency of 0.92822 ± 0.00006. The integrated efficiency is evaluated after having 
weighted the pt and rj distributions of the additional CMU muons to be equal to that of 
muons from heavy fiavor decays in the simulation. 

In the simulation, the probability that a CMUP pair passes the SVXII requirements 
described in Sec. Ill is 0.2206 ± 0.0047. This efficiency is measured in the data using muon 
pairs acquired with the J/ip trigger. We use CMUP muons withp^ > 3 GeV/c and \ri\ < 0.7. 
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FIG. 16: Invariant mass distribution of CMUP muons paired with other CMUP muons in the event 
with (a) or without (b) a LI CMU primitive. In order to derive the efficiency from the numbers of 
J/ip candidates in plots (a) and (b), the histogram in (a) has an entry for each CMUP leg. Solid 
lines represent the fits described in the text. 



The efficiency is evaluated in two steps. For each event, we first randomly choose a CMUP 
muon. After weighting the Zq and rj distributions of these muons to be equal to those of 
CMUP muons from simulated heavy fiavor decays, we derive the SVXII efficiency ei from 
the number of CMUP muons that pass or fail the SVXII requirements by fitting the dimuon 
invariant mass distribution with a straight line plus two Gaussian functions (see Fig. 17). 
For events in which the first randomly chosen muon passes the SVXII requirements, we 
derive the SVXII efficiency €2 from the numbers of second muons that pass or fail the 
SVXII requirements. After weighting Zq and f] distributions of the second muons to be 
equal to those of CMUP muons from simulated heavy fiavor decays, we fit again the dimuon 
invariant mass distributions with a straight line plus two Gaussian functions (see Fig. 18). 
The probability that a muon pair passes the SVXII requirements is ei -£2 = 0.2365 ± 0.0013. 
This measurement of the SVXII efficiency rests on the verified assumption that the event 
vertex ^-distribution is the same in the data and the simulation (see Fig. 19). 

In the simulation, the efficiencies of the LI and L2 triggers are 1 and 0.9976, respectively. 
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FIG. 17: Invariant mass distribution of CMUP muon pairs in which a first randomly chosen muon 
(a) passes or (b) fails the SVXII requirements. Solid lines represent the fits described in the text. 




FIG. 18: Invariant mass distribution of CMUP muon pairs in which a first randomly chosen muon 
track satisfies the SVXII requirements and the second muon track (a) passes or (b) fails the SVXII 
requirements. Solid lines represent the fits described in the text. 
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z (cm) 

FIG. 19: Distribution of the event vertex along the beam hne in the data (•) and in the heavy 
flavor simulation (histogram). 

By using muon pairs with CMU primitives acquired with the charm trigger, we measure the 
LI and L2 trigger efficiency to be 1 ±0.001 and 0.99943±0. 00045, respectively. The L3 trigger 
is not simulated. The L3 trigger efficiency is dominated by differences between the online and 
offiine reconstruction code efficiency The relative L3 efficiency for reconstructing a single 
muon identified by the offiine code has been measured to be 0.997 ±0.002 [23, 58]. However, 

Online algorithms are faster but less accurate than the offline reconstruction code. 
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TABLE VIII: Summary of efficiencies for reconstructing muon pairs from heavy flavor decays in 
the data and in the simulation. The last column indicates the corrections applied to the simulated 
efficiencies and used to derive .Acorr in Table VII. 



Source 


Data 


Simulation 


Corr. 


COT tracking 


(0.996 ± 0.006)2 


(0.998 ± 0.002)2 


1 ± 0.013 


CMUP acc. and eff. 


(0.5057 ± 0.0032)2 


(0.5235 ± 0.0022)2 


0.933 ± 0.014 


LI CMU primitives 


(0.92822 ± 0.0006)2 


(0.8489 ± 0.0026)2 


1.196 ±0.007 


Sili acc. and eff. 


0.2365 ± 0.0013 


0.2206 ± 0.0047 


1.072 ± 0.024 


LI eff. 


1± 0.001 


1 ± 0.001 


1 ± 0.0014 


L2 eff. 


0.99943 ± 0.00045 


0.9976 ± 0.001 


1.002 ± 0.001 


L3 eff. 


0.90 ± 0.01 


1 


0.90 ± 0.01 


Total 


0.0471 ± 0.001 


0.0435 ± 0.001 


1.084 ±0.035 



in a large fraction of the data, the L3 trigger has selected muons with the requirement that 
the distance between the track projection to the CMP chambers and CMP stub be Ar0 < 25 
cm, whereas the offline analysis requires Ar0 < 40 cm. Wc have measured the efficiency 
of this L3 cut by using J/ip candidates acquired with the J/ip trigger that has no Ar0 
requirement. After weighting the pr distribution of muons from J/ip candidates to model 
that of muons from b decays in the simulation, we measure the efficiency to be 0.948 ± 0.005 
for a single muon. The reconstruction efficiencies in the data and in the simulation are 
summarized in Table VIII. 



VII. DIMUON CROSS SECTION AND COMPARISON WITH PREVIOUS RE- 
SULTS 

We have selected pairs of muons, each with Pt > 3 GeV/c and \r}\ < 0.7, with invariant 
mass 5 < m^^ < 80 GeV/c2 and produced by double semileptonic decays of heavy flavors. 
The production cross section is given by 
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where N ^ BB ^ 52400 ± 2747 for bb production {N ^ CC ^ 19811 ± 2994 for cc 
production). The geometric and kinematic acceptance, .4corr, that includes trigger and 
tracking efficiencies measured with the data is hsted in Table VII. The integrated luminosity 
of the data sample is jC = 742 ± 44 pb~^. 

We derive 0-5^^ 5^^ = 1549 ± 133 pb. where the 8.6% error is the sum in quadrature of 
the 1.2% statistical error, the 2.9% systematic uncertainty due to the fit likelihood function, 
the 4.2% systematic uncertainty in the removal of the fake muon contribution, the 6% 
uncertainty of the luminosity, and the 3.2% uncertainty of the acceptance calculation. 

We also derive Uc^^^c^^ — 624 ± 104 pb. In this case, the statistical error is 6.4%, the 
uncertainty due to the fit likelihood function is ±8%, and the uncertainty in the removal of 
the fake muon contribution is 11.1%. 

We evaluate the exact NLO prediction of (75_>^ and (7c->^,c-+/i by complementing the 
MNR generator with the evtgen Monte Carlo program. We use the Peterson fragmentation 
function with e = 0.006 (0.06) for b (c) quarks and the measured fragmentation fractions 
listed in Table X. The NLO prediction is estimated using m — 4.75 (1.5) GeV/c^, the 
factorization and normalization scale Hr = fJ-F = -y/p|H-"m^, where m is the b (c) quark 
mass, and the MRST PDF fits [7] (we use the five fiavor scheme and A5 = 0.22 GeV/c^). 
In the following, we refer to it as the standard NLO calculation. We generate heavy flavor 
quarks with Pt > 2 GeV/c and \y\ < 1.3. 

The values of a^j^^ i^^ and ac-^/j^^c-^f^ predicted by the standard NLO calculation have a 
2% uncertainty, estimated by using different but reasonable procedures to sum the positive 
and negative weights, due to real and virtual soft gluon emission, returned by the MNR 
computation. The theoretical prediction also carries the uncertainty of the semileptonic 
branching fractions b ^ fi = 10.71 ± 0.22, b c ^ /i — 9.63 ± 0.44, and c ^ /i — 
9.69 ± 0.31% [60]. In the simulation, 79.4% of the muon pairs are due to 6 — > // decays, 
1.3% to 6 — > c — > /i decays and the rest to a mix of these decays. Therefore, the rate of 
predicted dimuon pairs due to bb (cc) production has a 3.7% (6.4%) uncertainty. For muon 
pairs selected with the same kinematic cuts of the data, the standard NLO prediction is 
Cb^t^fi^n = 1293 ± 55 pb. For cc production, the standard NLO prediction is (Jc^^/--^^ = 
230 ± 16 pb The ratio of the data to the standard NLO prediction with the above 

For comparison, the prediction of the herwig generator is af,^^ = 904 ± 33 pb and C7c^^,c^/a = 
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mentioned uncertainties {Rl — 1.20 ± 0.11 for bb production and Rl — 2.71 ± 0.49 for cc 
production) can be used to extract the value of cr^^ and acc in the data. In addition, as 
discussed in Sec. I, the theoretical prediction has a 15% uncertainty due to the choice of 
the heavy quark pole-mass PDF fits, and renormalization and factorization scales. After 
including the latter uncertainty, the ratio of the data to the standard NLO prediction is 
1.20 ± 0.21 for bb production and 2.71 ± 0.64 for cc production. 

The cc correlation measurement has no previous result to compare with. However, the 
CDF study in Ref. [25] has measured the ratio of dimuon pairs due to cc to those due to bb 
production. That study uses muon pairs selected in the same kinematic region as our mea- 
surement, and finds a ratio CC / BB = 0.15±0.02 whereas our fit to the impact parameter 
distributions yields CC / BB = 0.38 ± 0.07. In the simulation, this ratio is 0.17 (0.16) when 
using the HERWIG (mnr) generator. 

The extraction of from the dimuon production cross section and the comparison to 
other measurements is not a trivial issue. Muons with pt > 3 GeV/c and \r)\ < 0.7 are 
mostly contributed by b quarks with px > 6.5 GeV/c and \y\ < 1. However, there are 
tails contributed from b quarks with pt as small as 2 GeV/c and \y\ as large as 1.3. If 
these contributions are included, the resulting value of cif,^ is dominated by the production 
of b quarks with the smallest pr that, unfortunately, has a large statistical error because of 
the small kinematic acceptance. The measurement of a^i in Ref. [22] is based upon muon 
pairs selected with the same kinematic cuts as this study. That study does not report the 
value of a"{,^^5^^ but, in the assumption that these muon pairs are produced by b quarks 
with Pt > 6.5 GeV/c and \y\ < 1, quotes af,i{pT > 6.5 GeV/c, \y\ < 1) = 2.42 ± 0.45 /xb. 
It seems more appropriate to derive this cross section assuming that the ratio of data to 
theory is the same as that of the measured to predicted dimuon cross section. Using this 
method, the ratio Rl yields (Tiji{pt > 6.5 GeV/c, \y\ < 1) = 1324 ± 121 nb (the standard 
NLO prediction is 1103 ± 169 nb). In the simulation, only 75% of the muon pairs arise from 
b and b quarks with px > 6.5 GeV/c and \y\ < 1. The result of Ref. [22], rescaled by 75%, 

173 ± 11 pb. 

Following tradition, we vary the pole mass of b quarks by ±0.25 GeV/c^ and that of c quarks by 
±0.2 GeV/c^. 

The error is statistical. Systematic effects due to the fit likelihood functions or c-hadron lifetime were not 
investigated. 
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becomes (t,,i{pt > 6.5 GeV/c, \y\ < 1) = 1.80 ± 0.34 fjh. 

The op collaboration [24] has measured (T;,^^5_^^ = 1027 ± 260 pb using muon pairs 
with 6 < m^^ < 35 GeV/c^. That study selects muons with 4 < < 25 GeV/c, \ri\ < 0.8, 
and contained in a jet with transverse energy > 12 GeV. Reference [24] compares data 
to the exact NLO prediction that, evaluated with the hvqjet Monte Carlo program [61], 
is 357 pb. For this kinematical selection, except the request that muons are embedded in 
jets, the standard NLO prediction is (7b_»^ = 550 pb. When applying these kinematical 
cuts to our data, and before asking that muons are contained in jets with Ex > 12 GeV, we 
measure Uf,^^ = 658 ± 55 pb. 

Using the ratio Rl, the data yield cri,i{pT > 6 GeV/c, \y\ < 1) = 1618 ± 148 nb. The 
standard NLO prediction is af^{pT > 6 GcV/c, \y\ < 1) = 1348 ±209 nb ^\ For comparison, 
the HERWIG parton-level prediction is cri,i{pT > 6 GeV/c, \y\ < 1) = 1327 nb, and the cross 
section returned by the mc@nlo generator is a,,i{pT > 6 GeV/c, \y\ < 1) = 1704 nb, 27% 
larger than the mnr result 

The value of a^i has been extracted from the data using a fragmentation model based 
on the Peterson function. As previously noted, the mnr and HERWiG generators predict 
the same parton-lcvcl cross section Cf^ilpx > 6 GcV/c, \y\ < 1). The HERWIG generator 
models the 6-quark fragmentation differently, and this difference results in a prediction 
of (75_^^ which is 30% smaller than that of the mnr generator implemented with the 
Peterson fragmentation model. In the transverse momentum range of this study, the FONLL 
prediction for the single 6-quark cross section is fairly well reproduced by the nde calculation 
when using the Peterson fragmentation function with e = 0.002 [11, 62]. When using this 
fragmentation function, the exact NLO prediction becomes cr^^^^^^ = 1543 pb, which is 
20% higher than the standard exact NLO prediction. 

As argued in Ref. [63], the charmed quark production in e+e" data can be described at 
NLO accuracy using the Peterson fragmentation model with e = 0.02. In this case, the 

^"^ For charmed quarks, the standard NLO prediction is (JcciVT > 6 GeV/c, \y\ < 1) = 2133 ± 323 nb. 
We input the same 6-quark mass, scales and PDF fits used in the standard NLO calculation. 
The total hh cross section predicted by both mc@nlo and mnr generators is 56.6 fih and compares 
well with the result of the nde calculation (57.6 /xb). However, the inclusive single h cross section for 
Pt > 6 GeV/c and \y\ < 1 predicted by the nde and MNR programs are 5.5 and 5.6 fih, respectively, 
whereas mc@nlo generator predicts 11.8 /ib. In contrast, for both b and b quarks with pT > 25 GeV/c 
and \y\ < 1.2 the mc@nlo prediction is approximately 12% smaller than that of the MNR generator. 
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NLO prediction becomes ac-,fj,,c-*n — 383 pb, 66% larger than the standard NLO prediction. 
When using the smaller values of the e parameter, the ratio of data to theory becomes 
1.0 ± 0.2 for ai,^^ i^^ and 1.6 ± 0.4 for Gc-^^^c^^. 

VIII. CONCLUSIONS 

We have measured the production cross section of muon pairs from double semileptonic 
decays of h and b quarks produced at the Tevatron Fermilab coUider operating at ^/s — 
1.96 TeV. We select muons with pt > 3 GeV/c and \ri\ < 0.7. We select dimuons with 
5 < TO^^ < 80 GeV/c^ to reject the contribution of sequential decays of single b quarks 
and decays. The main sources of these muon pairs are semileptonic decays of b and 
c quarks, prompt decays of quarkonia, and Drell-Yan production. We determine the hb 
content of the data by fitting the impact parameter distribution of muon tracks with the 
templates expected for the various sources. Previous measurements of the bb correlations at 
the Tevatron yield contradictory results. The ratio of the data to exact NLO prediction is 
approximately 1.15 ± 0.21 when b quarks are selected via secondary vertex identifications, 
whereas this ratio is found to be significantly larger than two when identifying b quarks 
through their semileptonic decays. 

We measure a^^^ i^^ = 1549 ± 133 pb. The exact NLO prediction is evaluated using 
the MNR calculation complemented with the evtgen generator. In the calculation, we use 
rrib — 4.75 GeV/c^, the factorization and normalization scale fJ^R — fJ>F — \Jpt + "^6; 
the MRST PDF fits (we use the five flavor scheme and A5 = 0.22 GeV/c^). We use the 
Peterson fragmentation function with e = 0.006, and the PDG values for the fragmentation 
fractions. The NLO prediction is (Jb^^ ^^^ = 1293 ± 201 pb. The ratio of the data to the 
NLO prediction is 1.20 ± 0.21. 

From this measurement, we also derive (t^i{pt > 6 GeV/c, \y\ < 1) = 1618±148nb 
(the exact NLO prediction is 1348 ± 209 nb). The extraction of a^^i from the data depends 
on the choice of the fragmentation functions that connect a muon to the parent b quark. 
No fragmentation functions are available that match the accuracy of the NLO calculation. 
Reasonable changes of the fragmentation model indicate that the value of af^i extracted from 
the data has an additional uncertainty of approximately 25%. 
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APPENDIX A: SETTINGS OF THE HERWIG MONTE CARLO PROGRAM 

We generate generic 2 — > 2 hard scattering, process 1500, using version 6.5 of the HERWiG 
Monte Carlo program. In the generic hard parton scattering, bb and cc pairs are generated 
by HERWIG through processes of order (LO) such as gg bb (direct production). Pro- 
cesses of order are implemented in HERWiG through flavor excitation processes, such as 
gb gb, or gluon splitting, in which the process gg — > gg is followed by 51 — > bb. We 
generate final state partons with pr > 5 GeV/c^ and \y\ < 1.7. The hard scattering cross 
section is evaluated using the MRST fits to the parton distribution functions [7] . Hadrons 
with heavy flavor, produced by the HERWiG generator, are decayed with the evtgen Monte 
Carlo tuned by the BABAR collaboration [39]. We retain simulated events that contains 
a pair of muons, each of them with pt > 2.8 GeV/c^ and |?7| < 0.8 Wc find one good 
event in approximately 10^ generated events. These events are used to determine the kine- 
matical and detector acceptance as well as the impact parameter templates used to extract 
the heavy flavor composition of the data. Since different b hadrons, and especially different 
c hadrons, have quite different lifetimes and semileptonic branching fractions, it is impor- 

We also produced simulated samples requiring the presence of only one muon or no muons at all. 



43 



TABLE IX: Parameter settings used in our simulation are compared to the HERWIG default values. 



Parameter 


Default 


This Study 


QCDLAM 


0.180 


0.18 


RMASS(4) (c quark) 




1.50 


RMASS(5) (6 quark) 




4.75 


RMASS(13) 


0.75 


0.75 


CLMAX 


3.35 


3.75 


CLPOW 


2.00 


1.06 


LCLPW 




2.20 


DCLPW 




1.30 


PSPLT(l) 


1.00 


0.50 


PSPLT(2) 


1.00 


1.10 


CLSMR(l) 


0.00 


0.00 


CLSMR(2) 


0.00 


0.40 


PWT(3) 


1.00 


0.70 


PWT(7) 


1.00 


0.45 


SNGWT 


1.00 


1.00 


DECWT 


1.00 


1.00 


REPWT(O.l.O) 


1.00 


10.00 



tant that the generator models correctly the known fragmentation fractions and functions 
of b and c quarks. The HERWiG generator makes use of a large number of parameters that 
can be adjusted to this purpose. Unfortunately, their default setting [46] does not yield a 
satisfactory modeling of the heavy quark fragmentation that we have studied by comparing 
simulated decays (process 2160 of HERWIG) to e^e^ data. The available parameters do 
not allow to tune the ratio of baryon to mesons simultaneously for bottom and charmed fla- 
vors, nor to reproduce the measured ratio of vector to pseudoscalar resonances produced in 
the heavy quark hadronization. The first deficiency becomes a problem when the simulation 
of a QCD process, such as ours, is extremely time consuming. The second deficiency impacts 
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TABLE X: Fragmentation fractions in the tuned HERWIG simulation are compared to data. The 
fragmentation fractions of b quarks (first three rows) are defined according to the PDG notation [47] . 





Data 




HERWIG 


fu — fd 




(39.7±1.0)% 


39.6% 


fs 




(10.7±1.1)% 


11.2% 


/baryon 




(9.9±1.7)% 


9.6% 


f(c^D+) = f{c^D^) 




(16.4±2.3)% 


16.9% 


f{c^D*+) = fic^D*^) 




(22.8±2.5)% 


22.5% 


f{c^Ds + D*) 




(12.1±2.5)% 


11.5% 


f{c^ baryons) 




(9.5±4.0)% 


9.7% 



the evaluation of the kinematical efficiency and lifetime templates for c quarks because D* 
mesons mostly decay to mesons, the lifetime and semileptonic branching fractions of 
which differ by a factor of three from that of mesons. We have solved these issues by 
adding two additional parameters, analogous of CLPOW. In the HWUINC.F routine of the 
HERWIG program, the parameter CLPOW tunes the invariant mass distribution of cluster 
generated in the heavy quark hadronization. We use this parameter for b quarks only. For 
c quarks, the parameter CLPOW is replaced with two parameters, LCLPW and DCLPW, 
that separately control the yield of c-quark (mesons) and c-diquark (baryons) clusters, re- 
spectively. Table IX hsts the herwig parameter settings used in our simulation. Table X 
compares fragmentation fractions in the tuned simulation at the Z-pole to the data. The 
fragmentation fractions for b quarks are taken from Ref. [47]. For c quarks, the fragmenta- 
tion fractions are taken from Refs. [48-52]. The fragmentation functions are tuned in the 
HERWIG simulation by adjusting the parameters PSPLT and CLSMR to the values hsted in 
Table IX. Figures 20 and 21 compare some fragmentation functions predicted by the tuned 
HERWIG simulation to data. Simulated fragmentation functions are derived using decays 
generated with process 2160. Figure 20 compares the distribution of the fraction of energy 
of parent b quarks carried by all B hadrons resulting from the heavy quark fragmentation 
to OPAL data [53] that in turn are consistent with Aleph and SLD measurements [54, 55]. 
Figure 21 compares the fraction of momentum of the parent c quarks carried by D* mesons 
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FIG. 20: Distribution of ratio of the energy carried by all B hadrons to that of the parent b quarks. 
The data are OPAL measurements at the Z pole [53], while the herwig distribution is obtained 
with the parameter tuning listed in Table IX. 

to BELLE data [48] that in turn are consistent with the CLEO result [49] and the Aleph 
measurement at the Z-pole [56]. The data are fairly well modeled by the herwig generator 
with the parameter settings listed in Table IX. A similar agreement for the fragmentation 
functions, but not the fragmentation fractions, can be achieved using the tuning proposed 
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FIG. 21: Distribution of ratio of the momentum carried by D* mesons to that of the parent c 
quarks. The data are BELLE measurements [48], while the herwig distribution is obtained with 
the parameter tuning hsted in Table IX. 

by some of the HERWIG authors [57]. 
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APPENDIX B: RATE OF FAKE MUONS 

Muons reconstructed in the CMUP detector are divided in this study into real and fake 
muons. Real muons originate from semileptonic decays of hadrons with heavy flavor, the 
Drell-Yan process, and T decays. Fake muons include muons from tt or K decays and 
hadronic punchthroughs that mimic a muon signal. The probability that a tt or K track is 
misidcntificd as a muon is evaluated using — > Ktt decays reconstructed in data collected 
with the CHARM trigger. We select oppositely charged particles, each with pt > 3 GeV/c 
and I77I < 0.7, with \Szo\ < 0.5 cm. We require that each track is reconstructed in the 
microvertex detector with hits in at least four of the eight silicon layers. We evaluate the 
pair invariant mass for all pion-kaon mass assignments. The invariant mass is evaluated by 
constraining the two tracks to originate from a common point in the three-dimensional space 
(vertex constraint). We reject pairs if the probability of originating from a common vertex 
is smaller than 0.0002 or their invariant mass is outside the interval 1.77 — 1.97 GeV/c^. We 
also require that the displacement of the D'^-candidate vertex from the primary event vertex, 
projected onto the transverse momentum vector, be larger than 0.02 cm. To further 
reduce the combinatorial background, we also require the candidate to originate from 
a D*^ decay. We reconstruct D*^ decays by combining candidates with all additional 
COT tracks with a distance \6zq\ < 0.5 cm with respect to the vertex. Additional tracks 
are assumed to be pions and the D* invariant mass is evaluated by vertex constraining pion 
and D'^ candidates and rejecting combinations with probability smaller than 0.0002. The 
observed mD*± — rririo distribution is shown in Fig. 22 (a). We retain D° candidates with 
0.144 < mD*± — rriDO < 0.147 (their invariant mass distribution is plotted in Fig. 22 (b)). 
The fake muon probability is derived using the invariant mass spectrum of ttK decays 

in which one of the decay products is matched to a CMUP stub (see Fig. 23). We fit the 
data with two Gaussian functions to model the signal and a polynomial function to 
model the underlying background. The Gaussian functions model separately the right and 
wrong sign D° decays. In the fits to the data in Fig. 23, the width and peak of the first 
Gaussian function and the peak of the second one are constrained to the value returned 
by the best fit to the data in Fig. 22 (b) (peak at 1.865 GeV/c^ and a = 0.008 GeV/c^). 
Using the same method, we also evaluate the rate of fake CMUP muons that pass or fail the 
stricter < 9 selection cut described in Sec. VB (see Fig. 24). By using 361902 — > Kn 
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mp.±-mj30 (GeV/c^) m^o (GeV/c^) 



FIG. 22: Distributions of (a) m£,*± — m£)0 and (b) the invariant mass of candidates (the solid 
Une represents the fit described in the text). 



candidates we measure the fake muon probabilities listed in Table XI. The fake muon 
probabilities have been evaluated after weighting the transverse momentum distributions of 
kaons (pions) from I?° decays to model that of kaons (pions) produced by simulated 6-hadron 
decays (unweighted distributions are shown in Fig. 25). Since the fake muon probability is 
not a strong function of the transverse momentum (see Fig. 26), we ignore the statistical 
uncertainty of the simulated distributions because its effect is negligible compared to the 
10% uncertainty of the kaon and pion rates predicted by the simulation. 

TABLE XI: Probabilities, Pf and P/, that pions and kaons, respectively, mimic a CMUP signal 
for different selection criteria. 



CMUP selection Pf {%) P] V/o) 

standard 0.483 ± 0.003 0.243 ± 0.004 

< 9 0.347 ± 0.003 0.219 ± 0.003 

> 9 0.136 ± 0.001 0.025 ± 0.002 
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m^o (GeV/c^) m^o (GeV/c^) 

FIG. 23: Invariant mass distribution of candidates with (left) a kaon or (right) a pion leg 
identified as a CMUP muon. Solid lines represent the fits described in the text. 
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FIG. 24: Invariant mass distribution of D'^ candidates with (left) a kaon or (right) a pion leg 

identified as a CMUP muon. Top (bottom) plots require muons passing (failing) the < 9 cut. 
Solid lines represent the fits described in the text. 
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FIG. 25: Transverse momentum distributions of (left) kaons and (right) pions from decays 
are compared to those of tracks (•) arising from simulated 6-hadron decays. 
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FIG. 26: Fake muon probability as a function of the (left) kaon and (right) pion transverse mo- 
mentum. 
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